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Abstract 

The equation studied is u" + ^^^u' + eu u' + ku'^ = 0, with boundary 
conditions u{\) = 0, u{oo) = 1. This model equation has been studied by 
many authors since it was introduced in the 1950s by P. A. Lagerstrom. We use 
an elementary approach to show that there is an infinite series solution which 
is uniformly convergent on 1 < r < oo. The first few terms are easily derived, 
from which one quickly deduces the inner and outer asymptotic expansions, 
with no matching procedure or a priori assumptions about the nature of the 
expansion. We also give a short and elementary existence and uniqueness 
proof which covers all e > 0, /c > 0, and n > 1. 



1 Introduction 

The main problem is to investigate the asymptotics as e ^ of the boundary value 
problem 

Tl — 1 

u" + u' + euu' + ku'"^ = (1) 

r 

with 

= 0,u(oo) = 1. (2) 

We consider the cases k = and k = 1. Our interest in these problems, originally 
due to Lagerstrom in the 1950s [6], [7], was stimulated by two recent papers by 
Popovic and Szmolyan |10],[Tl], who adopt a geometric approach to the problem 
when k = 0, and there are many papers which use methods of matched asymptotics 
or multiple scales, with varying degrees of rigor. We will review some of this work 
below. The point of this paper is to give a completely rigorous and relatively short 
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answer to the problem without making any appeal either to geometric methods or to 
matched asymptotics. We can express the solution as an infinite series, uniformly 
convergent for all values of the independent variable. From this series we obtain the 
inner and outer asymptotic expansions with no a priori assumption about the nature 
of these expansions. An important, and as far as we know, original, feature is that 
there is no "matching" . 

Lagerstrom came up with these problems as models of viscous incompressible 
{k = 0) and compressible {k = 1) fiow, so much of his work centered on n = 2 or 
3, but he also discussed general n > 1 [8J. The infinite series we develop can be 
obtained for any real number n. What n controls is the rate of convergence of the 
series. 

For e = k = 0, there is an obvious distinction between n > 2 and n < 2. If 
n > 2, then the problem ([I]) — ([2]) has the unique solution 



so that the solution with e small is presumably some sort of perturbation of this. 
If n < 2 then there is no such solution. A consequence is that the convergence as 
£ — > is more subtle when n < 2 then when n > 2. Our analysis will show that 
there is little prospect of discussing the behavior for small e if n < 2, but fortunately 
we can handle all n > 2. Although it has been thought that finding the asymptotics 
when A; = 1 is considerably more difficult than when k = 0, [3], we will show that 
our technique covers each case with comparable effort. 

Our methods are not restricted to Lagerstrom's problems ([I]) — ((2]). In subsequent 
work (in preparation), we will show that there is a general method which can yield 
similar results for a class of singularly perturbed boundary value problems. 

We start in section [2] by showing that each of these problems has one and only 
one solution, for any n > 1 and any e > 0. This is based on a simple shooting 
argument plus a comparison principle. These results have been obtained before, 
but our proof is quite short. In the subsequent sections we develop the integral 
equation referred to above, and show how it leads with relative ease to the inner 
and outer expansions. These expansions go back to Lagerstrom and Kaplun, with 
rigorous justification of some of the features to be found in [1] or [11], for example. 
We find the exposition in Hinch's book f3j particularly clear (though nonrigorous). 
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and make that our point of comparison in checking that we get the same expansions 
as were found previouslylll 



2 Existence and uniqueness 

As far as we know, the first existence proof was by Hsiao [Ij, who only considered 
n = 1 and sufficiently small e > 0. Subsequently Tam gave what seems to be the 
first proof valid for all e > and k > 0, [Hj. Subsequent proofs by MacGillivray 
Cohen, Lagerstrom and Fokas [T], Hunter, Tajdari, and Boyer [5], each of which 
covers all e > 0, and by several other authors, e.g. [I2],|in], for restricted ranges 
of e, add to the variety of techniques which have been shown to work. Uniqueness 
is proved in [5J (for A; = 0) by use of a contraction mapping theorem, and in [Ij 
by essentially a comparison method. The goal of is not to give a short proof, 
but to illustrate the application of geometric perturbation theory to a much studied 
problem in matched asymptotic expansions. The proofs we give of existence and 
uniqueness are considerably shorter than the others we have seen. 

Theorem 1 There exists a unique solution to the problem ([21) — (III) /^'^ ^^?/ k > 0, 
e > 0, and n > 1. 

Proof. Like some others, starting with [T^, we prove existence using a shooting 
method, by considering the initial value problem 

Tl — 1 

u" + u' + euu' + A;m'2 = (4) 

r 

m(1) = 0, u'{1)=c, (5) 

for each c > 0. Since u' = implies that u" = and u is constant, any solution to 
this problem is positive and increasing. As was observed in 

u" + euu' < 0, 

and so from ([5]) , 

/ 1 2 

u + -eu < c. 

^We thank the referees for some very helpful comments. In particular, they called our attention 
to earlier proofs of the existence and uniqueness results, in some cases by techniques similar to ours. 
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In particular, since u' > 0, 

^ < , (6) 

so the solution exists, and satisfies this bound, on [l,oo). Therefore, limr^ooU (r) 
exists. Writing the equation in the form 

(r^-^w')' + {eu + ku') {r'^-\') = 0, (7) 

and integrating twice, gives 

r^-iy' (r) = ce-*^"W-^i'>W'^^ (8) 



u r 



c 



If M (2) < 1, then since u is increasing, implies that 
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C 



n(2)>p(c) = —^e-'-'ds. 

From this, and , we see that there are ci and C2, with < ci < C2, such that if 
c = Ci then u (cxd) < 1, while if c = C2, then u (oo) > u(2) > 1. Further, from ([8]) 
for any r > R > 2, 



u{r) = u{R) + c [ ^e-^"-^i° ""^*ds. 



If n > 2 the second term is bounded above by („_2)^i^n-2 ; while if 1 < n < 2 and 
i? > 2, it is bounded by c e~^^^~'^^^^'^^ds. Hence, this term tends to zero as i? ^ oo, 
uniformly for r > R, Ci < c < €2- Since u (R) is a continuous function of c, for any 
R, it follows that -u (00) is also continuous in c, and so there is a c with m (00) = 1, 
giving a solution to ([1]) — ([2]) . 

For uniqueness, suppose that there are two solutions of ([1]) — ([2]) , say ui and 
U2, with u[ (1) > u'2 (1) > 0. Then ui > U2 on some maximal interval, say {1,X) 
where X < 00. For the same initial conditions, if e = k = then direct integration 
shows that ui > U2 on (l,oo), and moreover, ui (00) > U2 (00). We then raise e 
and fc, looking for a pair (ei, /ci) such that ui (X) = U2 (X) for some X < 00, and 
if < £ < £1 or < A; < /ci, no such X exists. Hence, at {ei,ki) , Ui > U2 on 
[0, 00). If, at (ei, ki) , X < 00, then -Ui and U2 must be tangent at X, since Mi — U2 
has a minimum there, contradicting the uniqueness of initial value problems for ([1]). 
Hence, X = 00, and Ui > U2 on (1, 00). 
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Observe from ([7]) that if u'l (r) = U2 (r) for some r, then (r"^^ {u'l — u'2))' < 0, 
since ui > U2, so that there cannot be oscillations in u'l —u'2. Hence, ui (00) = U2 (00) 
implies that there is an R with u[ (R) = u'2 (R) and u'^ < u'2 on [R, 00). Integrating 
([7]) , and recalling that Ui{oo) = ^2 (c)o) , gives 

r-H<-U2)\^ = leR-' {ul-ul)\n 



-1 /»00 POO 

(n - 1) / (^2 _ _ ^ / ^n-l (^^/2 _ ^/2^ 

2 Jij Jr 



+ 2^(n 

The left hand side is zero, and all the terms on the right are positive, giving the 
necessary contradiction. ■ 

Remark 1 The existence theorem in IT^ has one added part. It is shown there that 
as 6 —>■ 0, the solution tends to a so-called "singular" solution obtained by taking 
a formal limit as e —>■ 0. See 110] for details. This limit result follows from our 
rigorous asymptotic expansions given below. 

Remark 2 There would seem to be no difficulty in extending the existence proof 
even to n < 1, but the uniqueness proof does use essentially the fact that n > 1. 

3 The infinite series (with /c = 0, n > 2) 

Starting again with ([T]), and u (1) = 0, we first consider the case k = 0, and obtain 



for some constant B. Since u (oo) = 1, ( ITOl) implies that u' (r) is exponentially 
small as r ^ oo. Hence we can rewrite (ITOl) as 



so that 



u — 



1 = C [ ^e-"*-"^-("-^)'^'rft. 

J oo ^ 

Setting er = p, et = r, and es = a, we obtain 

u{p)-l = r J_e-rg- ^^^^ 



oo ' 



where we use the arguments p and a to indicate that we mean the rescaled version 
of u. Here C is a constant satisfying 



-1 = C'£"'2/" ^e-"e-^-("-^)^"rfr. (12) 



T 



Since for each e there is a unique solution, this determines a unique C, dependent 
on £. 

We now consider the integral 

/oo 
{l-u{a))da, (13) 

which appears in the exponent in f|T2l) . This integral has been seen to converge for 
each £, but we need a bit more, namely, that it is bounded uniformly in e < r < oo 
as e ^ 0. To see this, we note that as a function of cr, u satisfies 

(Pu — 1 du ^ du Q 
da"^ a da da 

u = when a = e^u (oo) = 1. 

Denoting the unique solution by (cr) , we claim that if < ei < £2, then > 
for 62 < a < 00. If this is false, then ei and 62 can be chosen so that u^-^ (ctq) = 
(o"o) for some ctq > £2- But then, the problem 

d'^u n — 1 du du ^ 
da'^ a da da 
u (o-q) = Mei ((Xo) , u (00) = 1 

has two solutions, contradicting our earlier uniqueness proof. 

A consequence of this is that JJ^ {1 — u^-^ (cr)) da < (1 — Ue2 (c)) da, which 
implies that the integral in the exponent in ( fT2l) , including the minus sign in front, 
is bounded below independently of r > 5, and of e. We then see that the r-integral 
in (fT2l) approaches — 00 as e — >^ 0, and hence, that 

lim Ce"-2 = 0. 
Since f ^ {u — 1) da > 0, it follows from (ITT!) that if 
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then 



\u{p)-l\<Ce^-'En-i{p) 



(14) 



For purposes of future estimates, we make the obvious remark that 

( O (p2-") as p ^ if n > 2 
{p)=\ O (logp) as p ^ if n = 2 
[ O {p^-'^e-P) as p ^ oo. 

Hence if n > 2 there is a constant K such that 

K_i(p) <Kmin (p2-,pi-"e~^). 



(15) 



(16) 



The method now is to work from (ITTl) . As observed before, since u' (r) is expo- 
nentially small as r — s> oo, the integral term J^°° {u — 1) da converges. Hence, for 
given e > and po > 0, and any p > po. 



u{p)-l = Ce 



n-2 



r 



n-1 



- 1) (io- + 



— 1) do" 



(17) 



where the series in the integrand converges uniformly for < r < oo. 



In fact, we will need to use this series for all p > e. Thus we need to check its 
convergence in this interval. This follows from ffT^ and (fT5|) , which imply that for 
any p > e, if n > 2, then 



(u (s) — 1) ds 



< Ce 



n-2 



En-i is) ds 



and 



.n-2 



En-l{s)ds = I 



o (1) as e if n > 2 
as e ^ if 72 = 2 



Hence for > 2 and any C, the series in the integrand of (jT7j) converges uniformly 
on [e, oo). 
Now set 



/oo 
En-1 (s) ds. 



We note that, if n > 2, then $ as e 0, while if n = 2, then $ — > as C ^ 0. 
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We proceed to solve (ITTl) by iteration. Thus, the first approximation is, from 



u 



(p) - 1 = Cs-'' r -^e-^dr + O ($2) 

J oo ^ 



and we obtain the second approximation by substituting this back in (IT7|) . Repeating 
this, we reach 



u 



H iCe-'f Z 7^e-^ Ul En-.dafdr (19) 
- (Ce-'-f Z 7^e-^ L {/: 7^--' UL En-idt) ds} dadr + O m , 
as $ ^ 0. 

To obtain C, we need to be able to evaluate each of these terms for small p (in 
particular, for p = e), and this is a matter of integration by parts. Thus, for 
non-integral n, 

En-i [p) = I -rdr 



^n— 1 



1-n 



2-n 



T 



n-2 



dr 



P 



2-n 



-e-P + 



1 



n-2, 



. - . . ^, (20) 

2 — n 2 — n 

and this can be repeated to give En-i as a sum of terms of the form CkP^e~'' and 
En~p, until < n — p < 1. Then 

-E'n-p = / dr — / c/r 



rU-p 



^ 

r (p + 1 — n) 



-dr 



and we can then continue to integrate by parts as far as we like. (If n is an integer, 
we will reach ^-^dr, which introduces a logarithm.) 

Thus En-i (p) can be expressed as a sum of terms of the form CkP^e~P , and so 
obviously the same is true of E'^_^ , with e~^^ in place of e~^. Also, 



En-i {r) dr 



n-l 



da I dr 



T 



pEn- 



n-l 
'n-2 



da 



+ 



dr 



(21) 
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so that En-idr can be expressed as the same type of sum. Hence the second 
term in (fT9l) gives a sum of terms of the form Ek (2p) and the third and fourth terms 
a sum involving Ek (3p) . 

We now carry the process through in the most interesting cases, n = 2, 3. 



4 The case /c = 0, n = 2 

When n = 2 we are interested in 

r°° 1 

El (p) = / -e-^dr 
Jo ^ 



e log p + e ^ log rdr 



p 

poo pp 

= —e~''\ogp+ / log rdr — / log rdr 

Jq Jo 

= —e~^ log p — J — p (log p — 1) e^^ + O (p^ log p) , for small p, 

= - logp - 7 + p + 0(p2 log p). (22) 

(See, for example, 0, Chapter 1.) Also, for future purposes, using (l20i) we obtain 

E2 (p) = — - El (p) (23) 
P 

= i + logp + (7-l)-^p + 0(p2logp) asp^O. (24) 
p z 



Looking now at (]T9|) , with p = e, we see that as e ^ 0, 

Cloge ^ -1 

and 

Hence the series in ( !T9|) is in powers of jj^- 
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Also, we will work our approximations (in order to compare the results with those 
of Hinch in [3]) to order — stty; so that (for example) 

0- (r) b (r) ^ / _3 1 

u = + -4^^ + O log =^ - 

log(-) log (i) V ^ 

for any fixed value of r (p of order e). This, as we shall see, necessitates finding 
C = — \tt < IH drr H tttt + ( log"^ | - 

log i-e) I log (i) log' i-s) V V 

and requires use of all the terms in (IT^ . 

With this in mind, we look at the second term of (1191) . Thus from 0211) , 

POO 

/ Eidr = -pEi + e-P, (25) 
J p 

so that the second term is 

r -e-^ [tE^ - e-^) dT = C^[ f e'^ Eidr - T '—^drX 

J OO ^ L oo J oo ^ J 

= C'!^[-e~-E,] l'oo-2£^rfr} 
= C {-e-PE, (p) + 2E, (2p)) . (26) 
From (122|) . the second term is therefore 

C2(logp + 7-21og2p-27 + 0(p)) 

= (- log p - 7 - 2 log 2 + O (p)) (27) 

as p — > 0. 

In the third and fourth terms of f|T9|) we need only the leading terms, i.e. we can 
ignore the equivalent of 7 + 2 log 2 in fl27|) . Using fl25l) the third term becomes 

Ic"" f —{e-^ -TE,fdT=\c\\ogp + 0{l)) asp^O. (28) 
2 Joo 2 
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Finally, in the fourth term, the integrand in the r-integral is just the second term, 
(as a function of a), so that from (1251) , the fourth term is 



M 



{-e'^Ei (a) + 2Ei (2ct)} da 



dr. 



(29) 



It is seen from f l2^ that for any r < oo, Ei (a) da converges. Hence we can 
write the inner integral above in the form + , and it follows that 



/ {K + r{T)}dT 



M 



where i^T is a constant, r is bounded and r (r) = O (r logr) as r — > 0. 
follows that 

M = C^ {KEi (p) + O (1)) as p ^ 0. 
We can evaluate K using (l25ll and (l22l) : 



It further 



Ei(2a)rfa = - E,{u)du = -, 

oo 

e-'^E, (a) da = [- (e"'^ - l) E^] ^ - I (e"'^ - l) 



-da 



a 



= lim {-^1 (2a) + E^ (a)} = lim (log 2a - log a) = log 2. (30) 

(T— >0 (T— >0 

Hence, from fl29l) . the fourth term of f|T9|) is 

C {El (p) (log2 - 1) + O (1)} = -C^ {(log2 - 1) logp + O (1)} as p ^ 0. (31) 
Now setting p = e and using fl27j) , fl28l) . and fl3T|) . we obtain that 

-1 = -C(-log£-7 + 0(£)) + C'(-log£-7-21og2 + 0(£)) 
+ ^C' (loge + O (1)) - C {(log2 - 1) loge + O (1)} 



as £ — > 0. Hence, 
1 



log (9 



C 1 - 



7 



log 



-CM 1 - 



7 + 21og2 

log (M 



+CM--log2 +0 lo 



1 

e 
(32) 
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and 



where 



1 A B 

^-log(l)+W(i)+log^(i)+ llog^(i) 



Hence, 



-7 + A - 1 = 0, 
+ (7 + 21og2) + ^ -log2 = 0. 



A = 7 + 1 

5 = 72 + 27 + 1 -log 2. 



Thus, for fixed r, p of order e, we have, with A = log (i) , 

1 A , 7 + 1 , (7 + if -I -log ^. ^ , 
M - 1 = I ^ + + I (logr + loge + 7) 

+ ^ f 1 + ^i2±ll^ (_ logr - loge - 7 - 2 log 2) 



AH A 
1 /3 



^ A3 V 2 



log 2^ (logr + loge) + O (A~^) 



so that, after cancellation, 

log r 7 log r 



u 



+ ^ + 0(A-) 



A A2 



This is the "inner expansion". For the "outer expansion", i.e. fixed p, r of order 1, 
we use (fTOjl , truncated to second order, to get 



1 = -E. (P) + ^) + ^ (2i?i (2p) - e-^E, (p)) + O (A"^) 



These results are in accordance with those of Hinch and of others on this problem. 
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5 The case /c = 0, n = 3 

Here we are interested in (from (!22|) and (l23l) ) 

g-p 1 1 

E2 (p) = ^1 (p) = - + logp + (7 - 1) - -p + O (p2 logp) as p ^ 0. 

P P 2 ' ' 

Thus, the first term on the right of (IT^ evaluated at p = e is 

-C (1 + eloge + (7 - 1) + O (e^)) as e 0. 

The second term is 



{Ce) 



E2do I dr 



-Eo (t) 



E2 (a) (icr 



00 
p 



J 00 «^oo 

p 



(C^)^ -E2(p)/ E2ir)dr+ E'^dr). 



From fl23l) we see that 



/ E^rfr = h log^ p + O (logp) as p ^ 0, 

^00 P 



while from (l2Ti) 



£'2C?T = p-E'2 - -El = 1 + log p + 7 + O (p log p) 

' 00 

r E2dT = 1 logp + + O (log' p) . 

Joo P P 



In all, the second term is 



^2 I 1, 7 + 2 

(Ce) <( --logp 



P 



P 



+ O (log' p) I . 



It is readily verified that the third and fourth terms in (fT9|) give O jC^e^ log' p 
which is negligible. Thus, evaluating (fT9l) at p = £, we have 

-1 = -Ce Q + log £ + 7 - l"] + {Cef (-^ loge - + O {C^e^ log' e) , 
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so that 



C = 1 - 2£ loge - £ (27 + 1) + O {e^ log^ e) 



Then, for fixed r, p of order e, we have 



u — 1 = —e (1 — 2elog£ — e (27 + 1)) ( — + loge: + logr + 7 — 1 

\er 

+ (-— (log e + logr) - +0(e^\og^e) , 

\ er er J ^ ' 



u = 1 - - -eloge ( 1 - - ) (logr + J +5(1 M _ 



+ 0{eHog'e). 
For fixed p, r of order e"^, we again use ( fT9l) . to give 
M - 1 = (1 - 2£ log£ - £ (27 + 1)) (p) 

+ £2 (p) (p) _ _ y ^2^^| ^ Q ^^3) _ 



(33) 



Again, these results are in agreement with those of Hinch, and others, although (133|) 
gives one term further. 

Remark 3. It is of interest to consider what happens when n < 2, since, at least 
for n > 1, there still exists a unique solution. The equation (fT9l) is still valid at 
p = e, but since En_i (p) is no longer singular at p = for n < 2, (fT9|) with p = e 
becomes merely an implicit equation for This tells us that C ^ 0, since 

e^'"^ 00, but we no longer get an asymptotic expansion. In particular, it is no 
longer obvious that C is unique. Of course, we know this from Theorem [1] if ra > 1. 



6 The case k = 1 

We can in fact treat a generalization which causes no further difficulties, 

T) — 1 

u" + - — -u' + f (u) u'^ + euu' = 0, (34) 
r 
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with the same boundary conditions. As before, we will compare our results with 
those of Hinch in [3] . 

As remarked in the proof of Theorem 1, the solution will necessarily have m' > 
so that conditions on / [u) are necessary only for < u < 1. We require only that 
/ be continuous and positive in this interval. 

Then (1341) can be written as 

+ / (m) u + eu = 0, 



so that 



/r 
udt + A 



for some constant A, where 

F{u)= I f (s) ds. 

This becomes 



or, on integration. 



G{u)-G (1) = C f -^e'''-'^^^''-^^'''dt, 

where 

Giu) = r e^^^^dv. 



In order to keep the manipulations simple and effect comparisons, we will consider 
from here the Lagerstrom model, where / {u) = 1, F [u] = u, G (u) = e" — 1. Then, 
with er = p, et = r, we have 



e — e 



r^e-"e-^-("-i)'''^(ir, (35) 

J oo 

and writing 



u-1 = (e - e) 
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we get 



-j-n— 1 
oo ' 



As in section [3], we can integrate by parts, and since < ^^^^ < 1 in < m < 1, 
we will develop a convergent series as before. To get the first three terms (necessary 
to give Hinch's accuracy when n = 2), we have from ( l36ll that 



e — e 



oo 



(37) 

As before, since — e ^ exponentially fast as p ^ oo, the series in the 
integrand converges uniformly for large r, so that (1571) is valid for large p. But 
again we need to extend it down to p = e. From flS^ we have 

e" - e < Ce^-'E^., (p) 

and so the convergence proof is the same as that preceding (I19p . 

Before proceeding further with n = 2, we make a couple of remarks about the 
simpler case n> 2. Then, as we saw in subsection O only two terms are necessary 
to give the required accuracy, and then fl37|) gives 



q-n—l I / pU p 

oo ' k ^ oo 

and since appears in what is already the highest order term, we can replace it 
by its limit as « — >• 1, i.e. K Thus we get, to the required order, 

e--e = -Ce^-^E^_^-[ / ^ / (e" - e) rfarfr. 

V ^ / J oo ^ J oo 

(We will proceed more carefully for n = 2.) This, apart from the factor is the 
same equation as we dealt with in section [5] (with e" — e in place of m — 1), and the 
solution can be written down from there. (If we had a general function / in place 
of 1, we would get 

F{u) _ F{1) ^ ^ _ '^^""^ r ( r (^Fiu) _ Fil)\ ^ \ ^ ^ 
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Turning now to the case n = 2, and F [u) = u, we need three terms on the right 
of (EZD . Thus, 

4^ = i--L(e--e) + 0(e"-e)^ as « ^ 1. (38) 
e" — e e /e^ 



We follow the method used just before (11911 and obtain from (1371) that 



- ^ {Ce-'f £ ;^e-- £ | £ ( £ E^^.dt^ dsj dadr + O ($^) 

= + Fi + F2 + F3 + F4 + O (<I'^) , say. 

As before, if n > 2 then this is valid for any C as e — 0, uniformly in p > £, while if 
n = 2, it is valid as C — 0. 

For n = 2 we can continue to follow the argument in section 4. Thus, as p — > 0, 

Fi = -C^ (- log p - 7 - 2 log 2 + O (p)) 

e 

F3 = ^C' (logp + 0(l)) 

F, = -\c' [(log2-l)logp + 0(l)]. 

The term did not appear before. Only the highest order term is needed for our 
expansion and this is 



P T 
P 



Now 



r Efda = [tEI] ^ + 2 r—E,dT 
Jo Jo 

POO 

= 2 Eidr = 2log2, from 
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Thus, 



F2 = ^C'E^ (log2 + O (plog^ p)) = -ic^ ^iQg2) logp + O (1)) , 



and, evaluating at p = e, we have 
l-e = C(log£ + 7 + 0(£)) 

- ^C^ (log £ + 7 + 2 log 2 + O (£)) + log £ (-2 log 2 + 1 - 2 log 2 + 2) + O (C^) , 

.i,c3r3-4iog2.or^Vor^Vo^^^ 



2e^ V \log^/ \log^/ \log£ 

Hence if 

e - 1 A B _ /, _4 ^ 1 



then 

\2 



ie - IV 

-7 (e - 1) + A - ^ = 0. 

e 

A = (76 + e - 1) , 



S - ^7 + (7 + 2 log 2) - ^^l^Jl+i^ (e - 1)3 (3 - 4 log 2) = 0. 

e e ze^ 

We can of course calculate B, but in fact its value will be be irrelevant to the level 
of approximation that we take. 

Then, for fixed r (p of order e ), we have, with / = log (^) , 

e"-e^(e-l)|y + ^^^^^~^ + ^^^^3~ | (logs + logr + 7) 

^-^^-^F-^l (- log^ - log^ - 7 - 21og2) 

+ (3 - 4 log 2) (loge + logr) + O (r^) 

^l-e+^^^^^ + ^-^^^^^ + 0{r'). (39) 

t t 
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(Note that the definitions of A and B were such that u = at r = 1 up to 
and including order so that to that order there can be only terms in logr, not 
constant terms. We do not need the explicit value of B.) To obtain u, we have to 
invert, so that 

r e — 1 ^fe — 1) 1 

u = \og< 1 + —j—logr + — ^-\ogr + 0{l-^) L 

For fixed p, r of order e^^, we have 
e"-e = (l + ^^^^) E, (P) + ^^^ (2i?i (2p) - e-^E, (p)) +0 (r^) . 



Thus 



e - 1 ^ l + ^--e \ El jp) ^ (e - 1)' i2E, (2p) - e-'^E^ (p)) 



e \ I J I I 

(e-lfEUp) 



2e^ P '^(^")- 

Again, these results are consistent with those of Hinch and others, except that Hinch 
has an algebraic mistake which in pUj) replaces 7 + 1 — ^by •y — 1 + -. 



7 Final Remarks 

Starting with Lagerstrom, the terms involving log e in the inner expansions have been 
considered difficult to explain. They are often called "switchback" terms, because 
there is nothing obvious in the equation which indicates the need for such terms, and 
because, starting with an expansion in powers of e, one finds inconsistent results 
which are only resolved by adding terms of lower order, that is, powers of e loge. The 
recent approach to the problem by geometric perturbation theory explains this by 
reference to a "resonance phenomenon" , which is too complicated for us to describe 
here [IO],[n]. 

In our work, the necessity for such terms is seen already from the equation fllip 
and the resulting expansion ([Tj 



u{p)-l = 06''-^ 1^ -^e-^ ll- j {u - 1) da + ^ (^j {u-l)da^ } dr. 
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In the existence proof it was seen in ([9]) that C = O (1) as £ ^ 0. On the right 
of (|T71) the first term is simply —Ce"'^^En-i (p) , and the simple expansions given 
for El and E2 show immediately the need for the logarithmic terms. There is no 
"switchback", because the procedure does not start with any assumption about the 
nature of the expansion, and there is no need for a "matching". 

A number of authors have noted that the outer expansion is a uniformly valid 
asymptotic expansion on [l,oo), and therefore it "contains" the inner expansion 
[5], though this is more subtle when k = 1 [13]. Our twist on this is that both 
expansions are contained in the uniformly convergent series defined implicitly by 
(ITTl) . The simple derivation of this series, via the integral equation (fTT!) is new, as 
far as we know. 
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